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> ■ Abstract 



We analyze homogeneous anisotropic cosmology driven by the dilaton and the self-interacting 
"massive" antisymmetric tensor field which are indispensable bosonic degrees with the graviton in 



m 
o 

the NS-NS sector of string theories with D-branes. We found the attractor solutions for this system, 



which show the overall features of general solutions, and confirmed it through numerical analysis. 
The dilaton possesses the potential due to the presence of the D-brane and the curvature of extra 



dimensions. In the presence of the non-vanishing antisymmetric tensor field, the homogeneous 
universe expands anisotropically while the D-brane term dominates. The isotropy is recovered as 



the dilaton rolls down and the curvature term dominates. With the stabilizing potential for the 
dilaton, the isotropy can also be recovered. 
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I. INTRODUCTION 



Low energy effective theories derived from the NS-NS sector of string theories contain 
the gravity, g^ u , the dilation, $, and the antisymmetric tensor field, B^ v . The existence of 
the last degree of freedom leads to intriguing implications in string cosmology l|. In four 
spacetime dimensions, the massless antisymmetric tensor field is dual to the a pseudo-scalar 
(axion) field [2J, and the axion-dilaton system is known to develop an unobserved anisotropy 
in our Universe, which can be diluted away at late times only in a contracting universe [3|. 

Such a disastrous cosmological situation can be resolved also in string theory which 
contains another dynamical object called D-brane [4|]. The gauge invariance on the D-brane 
is maintained through the coupling of the gauge field strength to the antisymmetric tensor 
field flflEJ. Then, the effective action derived on the D-brane describes the antisymmetric 
tensor field as a massive and self-interacting 2-form field. Cosmological evolution of such a 
tensor field has been investigated in Ref. [7] assuming that the dilaton is fixed to a reasonable 
value. Although the time-dependent magnetic B field existing in the early universe develops 
an anisotropy in the universe, it was realized that the matter-like behavior of the B field (B- 
matter) ensures a dilution of the anisotropy at late times and thus the isotropy is recovered 
in reasonable cosmological scenarios In such sense the effect of antisymmetric tensor 
field on the D-brane is distinguished from that of field strength of the U(l) gauge field 

In this paper, we investigate the cosmological evolution of the B-matter-dilaton system in 
our Universe which is assumed to be imbedded in the D-brane. The usual string cosmology 
with the dilaton suffers from the notorious runaway problem, which is also troublesome in 
the D-brane universe. In our study, the dilaton obtains two exponential potential terms 
due to the curvature of extra dimensions A, and the D-brane tension (the mass term of the 
3-matter) mg. It is interesting to observe that the dilaton can be stabilized for negative A 
9I which, however, leads to a contracting universe due to the effective negative cosmological 
constant in our Universe. When A is positive, the B-matter dominance will be overturned 
by the A dominance as the dilaton runs away to the negative infinity. As a consequence of 
this, the initial anisotropy driven by the B-matter can also be diluted away at late times. 

For a realistic low energy effective theory, string theory must be endowed with a certain 
mechanism generating an appropriate vacuum expectation value for the dilaton. In such a 
situation, the dilaton is expected to be stabilized at some stage of the cosmological evolution 
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affecting the dynamics of the antisymmetric tensor field. Taking an example of the dilaton 
stabilization, we will also examine the cosmological evolution of the B-matter and the dilaton 
in which the essential features of Ref. [7| are reproduced. 

This paper is organized as follows. In Section [III we describe the low energy effective ac- 
tion of the D-brane universe and the corresponding field equations. Before our main discus- 
sion, Section [Till is devoted to presenting homogeneous solutions in flat spacetime illustrating 
some interesting features of the B-matter-dilaton system in a simple way. In Section IIV[ we 
find semi-analytic and numerical cosmological solutions to observe an intriguing interplay of 
the curvature A and the D-brane tension 777.5. In Section [V] we consider the evolution of the 
anisotropic universe with the dilaton stabilization which leads to a satisfactory cosmology 
of the D-brane universe. We conclude in Section IVT1 



II. STRING EFFECTIVE THEORY ON THE D-BRANE 

The main idea of the D-brane world is that we reside on a Dp-brane imbedded in 10 (or 
11) dimensional spacetime with extra-dimensions compactified. The bosonic NS-NS sector 
of the D-brane world consists of the U(l) gauge field living on the Dp-brane and the 
bulk degrees including the graviton g^ u , the dilaton $, and the antisymmetric tensor field of 
rank-two B „. In the presence of the brane, the gauge inwia.ce of B„ is restored through 
its coupling to a U(l) gauge field and the gauge invariant field strength is [4J 

= B^ u + 2-Ka'F^, where = d^A v - d u A^. (1) 

Even though we are assumed to live on the D-brane, we adopt the conventional com- 

pactification in the sense that the extra dimensions are compact and stabilized, and thus 

static. In the presence of the D-brane, this would in general require an additional setup 

like additional branes and fluxes along the extra dimensions 10], and the working of non- 
j — I 

perturbative effects [11] . The extra dimensions are warped due to branes and fluxes. The 
warping of extra dimensions gives a warp factor in the definition of the four- dimensional 
Planck scale [12 ]. and also induces the potential for the dilaton [lOj as described below. If 
the fluxes in the extra dimensions significantly affect the compactification, the existence of 
the fluxes on the D-brane seems also natural and it is intriguing to tackle their effect in 
the early universe. However, simultaneous consideration of the fluxes along both the extra 
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dimensions and the D-brane lets the computation almost intractable. Therefore, the con- 
ventional compactification is an appropriate setup at the present stage, and then we take 
the following four- dimensional effective action of the bosonic sector in the string frame [7| 



5s 



1 



d 4 x 



-2$ 



R-2A + 4V^V^ - ^H pvp H^ p 



-m B e 



(2) 



where the tilde denotes the string frame quantity, H pvp = d[ p B up ] and B*^ u = \^f—gt pV apB afi 



2n\T p where T p is the effective 



with eoi23 = 1- The parameter mg is defined by m 
brane tension. Note that we omitted the Chern-Simons like term assuming the trivial R-R 
background. If we assume that the six extra-dimensions are compactified with a common 

1 r) p — 3 15— p 

radius R c , one finds m B = ^(gl/Air) ^ (R C M P )~*~ M P where g s = e* and M P = 2.4 x 
10 18 GeV is the four- dimensional Planck mass. The qualitative features of our results do not 
depend on specific values of p. Thus, we take p = 3 for simplicity. The A term comes from 
the scalar curvature or the condensate (H 2 ) of extra dimensions integrated over the whole 
extra dimensions. 

The action, and thus also field equations, can be written in a more familiar form in the 
Einstein metric, which is defined by 



9fiu 



-2<E>~ 



(3) 



We will work in this metric from now on. In terms of Einstein metric, the action becomes 



S 



2ni 



L — / <t'l\'~Il 



R - 2Ae 2 * - 2(V$) 2 - ^f~^H 2 



m|e 34 \/l + -e-^B 2 - — e" 8 * (B*B) 
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The field equations derived from the action (T4]) are 



B, 



/IV 
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\e~^B; u (BB*) 



1 + \e~^B 2 - ^e- 8 * (BB 



Ct - k 2 T 



(4) 



(5) 

(6) 
(7) 
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where the dilaton potential is 



2Ae 



2* 



4$ 
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H 2 



+ m|e 3 * 



1 + 2 6 



_4$ 



B 2 - 



-8$ 



16 



and the energy-momentum tensor is given by 



(B*Bf 



(8) 



-g^Ae 2 " + 2V M $V„$ - ^(W)' 



H e 
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.4$ 



1 2 w -^-5^ e " 4 * g2 + e " 8 * g ^ g . 

^ /ill _ — 4<J> K>9 1 __Sff) /*7k>*\2 



A 



l + |e- 



23 2 - — e" 

A 16 C 



(9) 



In the subsequent sections, we examine the equations of motion ©-([I]) and find homo- 
geneous solutions. Cosmological implication of the D-brane universe is of our main interest, 
including stabilization of the dilaton and condensation of the antisymmetric tensor field. 



III. HOMOGENEOUS SOLUTIONS IN FLAT SPACETIME 

For the sake of the simplicity, let us first look into the dynamics of our system in flat 
spacetime. With g^ u = rj^, the Einstein equations © become simple and gauge invariance 
of the field variables E l = B^ and B* = e r,fc£>jfc/2 which directly appear in the expression 
of energy-momentum Q allows a homogeneous configuration 

$ = E = E(t), B = B(t). (10) 

For homogeneous configuration (TTUj) . Oi-components (electric components) of the equation 
of the antisymmetric tensor field ([SD reduce to an algebraic equation 

E + e -«B(E.B) 

./I _ e -4* ( E 2 _ B 2) _ e -s<w E . B )2 



B - e - 4 *E(E ■ B) 



ij-components (magnetic components) give 

— 4B$ = —p 3 ® - = = ~ 7 — , (12) 

_ e - 4 *(E 2 - B 2 ) - e" 8 *(E • B) 2 

where the overdot denotes differentiation with respect to the rescaled dimensionless time 

t = m B t. (13) 



These equations (|TT]) - (|T2|) tell us that the magnetic components B are dynamical but the 
electric components are determined by the constraint equation f[TT|) . Dynamics of the dilaton 
is governed by §6§ which becomes 

* + ie--B 2 = -Ac" - V 3-e-^(E?-B') + e-»(E.B)« 
4 ' 4 ^/l - e -4*(E 2 - B 2 ) - e-^(E ■ B) 2 ' 

where the dimensionless parameter A is defined by 

A = A (15) 

The constraint equation (TXTT) forces its numerator to vanish except a trivial solution 
$ = — oo, however vanishing numerator allows only a vanishing electric field solution E = 0. 
Then the equations for the magnetic components (fT21) and the dilaton f[T4"j) reduce to 

B - 4B$ = -e 3 * - B =, (16) 
Vl + e-^B 2 V ; 

$ + V 4$ B 2 = -Ae 2$ - ^ d + e = . (17) 
2 4 Vl + e" 4 *B 2 

In the following, we will find solutions of (fl~6l) and (IT71) to see how the dilaton and the 

magnetic component behave. 



A. Dilaton potential with B = 

When the magnetic components of the antisymmetric tensor field B vanish, we have only 
the dilaton equation from (fi~7|) 

• = -»• (18) 

where the dilaton potential U(<&) is given by the sum of the curvature term and the brane 
tension term 

U(9) = ^ + \e™ (19) 

as shown in Figure [TJ 

For A > 0, the dilaton potential U(<&) monotonically increases, so the dilaton finally rolls 
to the minimum $(i) — > — oo as t — > oo. For A < 0, the dilaton potential £/($) has the 
minimum value i7 min = 8A 3 /27 at $ min = ln(— 4A/3). If E = t/ min , the dilaton is stuck at 
$ = $ m in- When U m i n < E < 0, the dilaton oscillates around the minimum. Therefore, the 
dilaton is stabilized with mass m$ = |A 3 / 2 m#, and due to exponential potential terms, its 
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FIG. 1: The graph of U(§). For A > 0, U(<&) monotonically increases (dashed line). For A < 0, 
U(<&) has a minimum at $ = ln(— 4A/3) (solid line). 

mass can be much different from the vacuum expectation value of the dilaton rriB^mm- As A 
increases, the dilaton mass rapidly grows. When E > 0, the dilaton again rolls to negative 
infinity. It means that the dilaton is easily destabilized by its fluctuations, if the dilaton 
mass becomes too small. Though the negative vacuum energy U m \ n < does not affect the 
dynamics of the dilaton and the antisymmetric tensor field in flat spacetime, its coupling to 
gravity leads to singular cosmological evolution as we shall see in Section IIV1 Therefore, for 
the proper description of cosmology, some additional term is needed to adjust the value of 
U min to vanish as we shall discuss in Section [Vj 

The runaway or stabilized behavior of the dilaton remains essentially unchanged even 
with B 7^ although the dependence on the initial conditions alters as we will see in the 
next subsection. 

B. Magnetic solution 

Without loss of generality, let us assume that the direction of the magnetic component 
of antisymmetric tensor field is fixed as B(t) = B(t)k. Although the equations of motion 
( TT6]) — (jTTT) are still complicated even under this assumption, the presence of the magnetic 
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components do not alter the story of the dilaton much. 

First, let us consider a simple case of zero brane tension (rag = 0), which is analytically 
tractable. Then in (jl6p - (ll7j) . only the A-dependent term survives on the right-hand side. 
(For this case, we will use the original A and time t since tub = 0.) The solutions to these 
equations are 

i T/1 _ \ 2 1 

+ \nc u (20) 



$(T) = T - - In 



and 



where & a , b 2 , and C\ are integration constants, and the time has been rescaled to T = Ci(t+c 2 ) 
by absorbing another constant Ci- The solutions are still remaining valid under T — > — T, 
i.e., under the change of the sign of c\. Therefore, we can fix the time T to flow to the 
positive direction and c\ > 0. The solutions approach their asymptotic configurations at 
large T, 



$(T) = -T + ln(4 Cl ), (22) 
±B(T) = signC&O^- + 6 2 = 5 st . (23) 

We observe that the dilaton rolls linearly in time to negative infinity, and the antisymmetric 
tensor field condensates to a constant ±-B st which depends on A. There are two topologically 
distinct condensation processes depending on the signature of b\ which is to be set by the 
initial conditions; the change in A contributes to B st oppositely. 

With the brane tension term turned on, we numerically study due to the complexity of 
field equations. The solutions are classified into two classes as in the pure dilaton case in 
Section UlIAl depending on the signature of A. 

For A > 0, the dilaton approaches negative infinity irrespective of initial conditions, $o 
and B . The condensed value of the magnetic component deviates from the initial value in 
the beginning, and then stabilizes to a constant as shown in Figure [2j We observed from the 
numerical results that there exist two branches of the magnetic condensation depending on 
the initial conditions as discussed above, but we show only one branch in the figure. Like 
in the mg = case, once the magnetic component is condensed, the condensation survives 
as a constant value. It is different from the case without the dilaton, in which the magnetic 
component permanently oscillates in flat spacetime 
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FIG. 2: and B in the Einstein frame for A = (solid lines), A = 1/2 (dashed lines), and A = — 1/2 
(dotted lines). The initial conditions are <&o = 1 an d -Bo = 0.1. The configurations are very similar 
regardless of A. 



B(t) 




400 



FIG. 3: <3? and B in the Einstein frame for A = — 1/2, B = 0.1, and $ 
stabilization is observed as well as the periodically waving oscillation of B. 



-0.2. The dilaton 



When A is negative, there are two classes of solutions, depending on the initial value of 
the dilaton $ ( or equivalently the energy density) for a fixed initial magnetic condensation 
Bq. For large $0 or Bo, the configurations are almost the same as those for nonnegative 
A (see Figure [2]). For small enough $0 and Bq, the dilaton shows the oscillating behavior 
as expected from the dilaton potential in Figure [U According to the dilaton oscillation, 



9 



the magnetic component is also oscillating around a condensed value (see Figure [3]). If we 
involve the radiation of perturbative modes of the dilaton and the antisymmetric tensor 
field, which is missing in this classical configuration, the classical solution will damp to a 
stabilized value of the dilaton and the final condensed value of the magnetic component. 



IV. COSMOLOGICAL HOMOGENEOUS SOLUTIONS 

In this section we study cosmological solutions of the D-brane universe in the presence of 
30th the dilaton and the antisymmetric tensor field. As we have done in the previous paper 
3], we assume spatially homogeneous configurations for the antisymmetric tensor field and 
the dilaton, and look for the time evolution of these fields and the expansion of the universe. 

The non-vanishing homogeneous antisymmetric tensor field, in general, implies the 
anisotropic universe. To consider the simplest form of anisotropic cosmology, we take 
only a single magnetic component of B^ to be nonzero, namely Bi 2 (t) = B(t) and 
Boi(t) = B 23 (t) = B 3 i(t) = 0. Then the metric consistent with this choice of field con- 
figuration is of Bianchi type I 

ds 2 = -dt 2 + ai (t) 2 (dx 1 ) 2 + a 2 (t) 2 (dx 2 ) 2 + a 3 (t) 2 (dx 3 ) 2 . (24) 

Then, in the Einstein frame, the field equation for B (J5]) reads 

d 2 B ( 1 dai 1 da 2 , 1 da 3 A d&\ dB m 2 B e^B 



dt 2 V «i dt a 2 dt a 3 dt dt J dt ^/l + B 2 /e^a 2 a 2 ' 



the dilaton-field equation ([6]) is 



d 2 $ / 1 dai 1 da 2 1 da 3 \ c?$ 1 

2q b - -Qb - Qb~ Qa, (26) 



dt 2 \ai dt a 2 dt a 3 dt J dt 2 

and Einstein equations ([7j) are 

1 da\ 1 da 2 1 da 2 1 da 3 1 da 3 1 dai 

+ —-JT — -JT + — -7T — -7T = 0* + QB + 26 + QA, (27) 



ai dt a 2 dt a 2 dt a 3 dt a 3 dt a\ dt 
1 d 2 a 2 1 d 2 a 3 1 da 2 1 da 3 

jl-f" H T^T H T+ TT = ~Q® + Qb + Qb + Qa, (28) 

a 2 dt z a 3 dt z a 2 dt a 3 dt 

1 d 2 a 3 1 d 2 a\ 1 da 3 1 da\ 

j7T H H T* 17 = + Qb + Qb + Qa, (29) 

a 3 dt z ai dt z a 3 dt a\ dt 

1 d 2 a\ 1 d 2 a 2 1 dai 1 <ia2 

+ — ~mT + — -JT — -JT = -69- Qb + Qb + Qa, (30) 



ai (it 2 a 2 dt 2 a\ dt a 2 dt 
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where 



(d*\ 2 e~ 4$ (dB\ 2 2$ 

1 2 3* A e" 4 *B 2 V /2 „ 1 2 3* A e - 4 *B^" 1/2 



(31) 
(32) 



12 / \ 1 2 

It is more convenient to employ the dimensionless time variable t of ffl3|) and to introduce 
the variables a* and b defined by 



oti = In a. 



aia 2 



Then the above equations (l25l) -( ]30il are rewritten as 



b + (6li + "2 + d 3 )5 + 



2$ + 2(-di - d 2 + d 3 - 2$)$ 



+di + d 2 + (di + d 2 )d 3 



3* 



$ + (di + d 2 + d 3 )$ = -2p B - p A - -p& - Pb, 



did 2 + d 2 d 3 + d 3 di 
di + di(di + d 2 + d 3 ) 
d 2 + d 2 (di + d 2 + d 3 ) 
d 3 + d 3 (di + d 2 + d 3 ) 



p$ + Pb + Pa + Pb, 

Pa + Pfe, 

Pa + Pb, 

2pB + Pa + Pb, 



where, by dividing m 2 B from each g in Eqs. f l3Tj) - fl32l . we have 



l r- • i 2 

p$ = $ 2 , p B = - |& + (d x + d 2 + 2$)6j , p A = Ae 2 *, 



Pb 



/2 



(33) 



(34) 

(35) 

(36) 
(37) 
(38) 
(39) 



(40) 



2 v / , ru 2 

where A is defined in ( fl5|) . In the subsequent subsections, we look for the solutions to the 
above field equations for various cases beginning with some simple solutions. 



A. Attractor solutions 



The dilaton in the system of equations (!34l) - (l39|) has the exponential potential up to the 
correction due to the antisymmetric tensor field. It is well-known that the scalar field with 



11 



the exponential potential possesses the scaling solution in which the energy density of the 
scalar field mimics the background fluid energy density [13I . Q]. This scaling solution is 
also an attractor, so that the late time behavior of the solutions are universal irrespective 
of initial conditions. This is an attractive property of the exponential potential. For the 
potential of the form V A ($) = Voe^*, there is an attractor solution 



2 

$ = --ln 



(3V Q 1/2 t 



%12-/3 2 ) 



2^ 2 



a= (-) hit, (41) 



for < (3 < v 12 . The scale factor obeys the power-law time-dependence, implying that the 
rolling of $ constitutes the matter having an equation of state p = wp where w = /3 2 /6 — 1 
varies from —1 to +1 for the aforementioned range of (3. 

We found this type of particular solutions of the Eqs. (!34l) - (l39l) . which can be found when 
we have a single exponential term in the potential, that is, for the case of A = and for the 
case of m B = 0. For both cases we start from an ansatz of the form 

$(t) = 7$ hit + $ , Oil = 7i hit, a 3 = 73 hit, b = constant, (42) 

where we suppressed constant terms in cti and 0:3 which correspond to simple rescaling of 
coordinates. For the case of A = 0, we obtain two distinguished solutions 

= ~ hit + In jj, ai{t) = a 3 (t) = - hit, b = 0, (43) 



and 



2 ~ 2 /8 ■ 5 1//4 \ ~10 3 1 

$(t) = __i n f + _hW__j ) ai(t) = -\nt, a 3 (t) = -\nt, b = ±-. (44) 



The first solution is nothing but the solution ff4T|) with (3 = 3. The second solution has 
the non- vanishing antisymmetric tensor field. For the case of = 0, we introduce a new 
rescaled dimensionless time variable i = A 1//2 t instead of t of (jT3]) and then get the continuous 
set of solutions from Eqs. (12"5]) -(13"0" |) 

= — lnt + - In 2, ai(t) = a 3 {i) = liai, b = arbitrary constant. (45) 

$(t) and a(t) are same as those in Eq. ( HIT) with (3 = 2, while we have the non-vanishing B 
field condensate. 

These solutions are the solutions to the Eqs. (!34]) - (l39|) for the specific initial conditions. 
However, the importance of these solutions, as noted in the first paragraph, arises from the 
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fact that they are attractors, which means that after enough time the solutions with different 
initial conditions approach these solutions. We will confirm this through numerical analysis 
in the next subsection. 

The solution (|43[) applies for the brane tension dominated case where the dilaton potential 
is approximated by m B U(<&) = ^m B e 3 ^ from Eq. (1191) . The evolution of the dilaton under 
this potential produces matter with the equation of state p = \p. 

Once the antisymmetric tensor field is turned on, the anisotropy appears as in the solution 
f f44|) . The measure of anisotropy is 

^ = -- (46) 

This result is contrasted with that in Ref. [7f where the dilaton is assumed to be stabilized. 
The rolling of the dilaton makes the difference. It affects the dynamics of B field in such a 
way that b remains constant at b = ±1/2 instead of oscillating about the potential minimum 
6 = and the anisotropy is maintained. 

Let us turn to the third solution (H5l) . It is relevant when the dilaton potential arising 
from the curvature of extra dimensions, V(<&) = Ae 2 *, dominates over other contributions. 
In our scheme this happens as the dilaton rolls down the potential. When b vanishes, the 
transition point at which Ae 2 * starts to dominate over |m|e 3lf is at $j = ln(2A). Thus, this 
solution describes the late time behavior of all the solutions with various initial conditions 
when A is positive. It is very intriguing since we achieve the isotropic universe in the end. 
The rolling of dilaton makes the brane tension term which causes the anisotropy when we 
have non-vanishing B field less important than the extra dimension curvature term which 
recovers the isotropy. The rolling of the dilaton under the potential Ae 2 * now forms the 
matter having the equation of state p = — |p, thus giving marginal inflation. 

B. Numerical analysis 

1. Initial conditions 

We have the second order differential equations for four variables $(£), Bit), = a2(t), 
a 3 (t). Thus we need eight initial values <3> , &o, B , B , a w , a Wl a 30 , «3o to specify the 
solution. Among these, ctio can always be set to zero by coordinate rescaling. d i0 must obey 
the constraint equation (1361) . but this does not fix the ratio dio/a3o- We choose the isotropic 
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universe as a natural initial condition which leads to dio = a 30- 

Since the dilaton potential is composed of exponential terms, the shift of the dilaton field 
by a constant can be traded for the redefinition of mass scale. We use this property to 
take the initial value of the dilaton to be zero without loss of generality. In our numerical 
analysis, we take the dimensionless time variable as t = fht where rh = m^ei* and use the 
variable $ = $ — $o with its initial value $o — 0. This means that the proper time scale for 
the cosmological evolution is not m^ 1 , but frT x . For convenience's sake, we take the initial 
time as % = 1. The other mass scale A is also affected by this shift and we can treat it by 
replacing the parameter A with A = Ae - * . 

Now we need three initial values $o> B Q , an d B , to fix the functional form of the solution. 
The initial values bo and bo are related to B and B by 60 = B and 60 = B — 2(di + $o)-Bo- 

2. Numerical solutions for B = 

If we assume B^ = 0, we can study the role of the dilaton interacting with gravity more 
clearly as was also true in the subsection IIII Al We keep the brane term with the tension 
tub- The spacetime is now isotropic. The parameter A in the string-frame action, which can 
be interpreted as the curvature of the compact internal manifold, plays a very important 
role, and the solutions are topologically different depending on the signature of A. 

When A = 0, we have an attractor solution ( 1431) . Since we consider the dilaton only, 
the only relevant initial condition is $o- The initial condition for the attractor solution is 
$0 — — a/3/2 in the setup described in the previous subsection. For a few other values of 
$0) we plot the numerical solutions in Figure HI Regardless of the initial value <3>o> all the 
solutions approach universally the attractor solution (up to rescaling for the scale factor) 
after some time. 

When A > 0, the evolution is divided into two stages. When $ is large so that p b = |e 3 * 
is much larger than pa = Ae 2 *, the solution approaches the attractor fj4*3l) . As $ rolls down 
to around $ t = ln(2A), p\ becomes larger than p b . After that point the solution approaches 
the attractor fT4"5l) with 6 = 0. In Figure [5], numerical solutions show this two stage evolution 
explicitly. 

Since A measures the curvature of the extra dimensions, its value can be negative. When 
A < 0, the dilaton potential has the global minimum at $ m i n = ln(— 4A/3) with negative 
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FIG. 4: Numerical solutions for B = and A = with initial value $o = — \/3/2 (solid curve, the 
attractor), (dotted curve), —3 (dashed curve) and 1 (dash-dotted curve), respectively. For the 
evolution of the scale factor, we plot ta(t) for it becomes constant 7 when the scale factor behaves 
as a oc i 7 . 

cosmological constant U min as discussed in the subsection IIII Al The disastrous consequence 
of this negative energy minimum is that the scale factor collapses in the end. This cannot be 
avoided because the dilaton rolls down toward the minimum and stays where the potential is 
negative while the kinetic energy is diluted by the expansion, thus the total energy becomes 
negative and it derives the universe to collapse. It is different from the case of flat spacetime 
in the subsection IIII A[ where the dilaton experiences either the permanent oscillation around 
$ m i n for negative initial energy density or the attractor solution for nonnegative initial energy 
density. 

3. Numerical solutions for B 7= 

Now, we turn on the antisymmetric tensor field along the x 3 direction, B12 = B 7^ 0. The 
spacetime becomes anisotropic, ai(t) = a 2 (t) 7^ 0,3 (t), in general. The solutions are classified 
again by the signature of A. 

For A = 0, we have an attractor (1441) under the initial conditions of $0 — —7/4 • 5 1 / 4 , 
60 = ±n> an d ^0 = 0. In Figure E], we plotted numerical solutions for a few different initial 
conditions in addition to the attractor. It is confirmed again that all the solutions approach 
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FIG. 5: Numerical solutions for B = and A > 0. When we set the parameter A = 10 , the curves 
stand for the solutions with the initial value $o = —^3/2 (solid curve) and (dotted curve), —3 
(dashed curve) and 1 (dash-dotted curve), respectively. 

the attractor as the time goes on. The final value of b is either +| or — | depending on 
the initial conditions. In the lower- right panel in Figure [6] we provide the evolution of each 
component of energy density. The kinetic energy of the dilaton p$ catches up the potential 
energy pb and the ratio of them becomes constant. This is a characteristic feature of the 
scaling solution [14]. The kinetic energy of B field is kept much smaller than both of them, 
but the anisotropy is still maintained due to the difference between pb and pV 

For A > 0, the evolution is again divided into two stages as in the B = case. In 
the first stage where pb is dominant, the solution approaches an attractor f|44l) . In the 
second stage where p\ is dominant, it approaches another attractor ([45]) . Thus the universe 
recovers the isotropy. The final value of b is a certain constant which is determined by initial 
conditions and can differ from ±|. Numerical solutions for a few initial conditions are shown 
in Figure [7] The kinetic energy of the dilaton p$ now catches up the potential energy pb in 
the first stage and p\ in the second stage. The ratios, p$/pb and p$/pa, approach constants 
in each stage. The kinetic energy of B field is kept much smaller as in A = case. 

For A < 0, the solution becomes singular as in the B = case. Here we skip the 
description of such singular solutions which are not suitable for the evolution of our Universe. 
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FIG. 6: Numerical solutions for B ^ and A = 0. The initial conditions are as follow. Solid line: 
$ = -7/4 • 5 1 / 4 , b = 1/2, B = 0; Dotted line: 6 = 0, b = 1, B = 0; Dashed line: 6 = 0, 
6o = 10, Bq = 0; Dash-dotted line: 60 = —^3/2, bo = 1, Bq = 0. The lower-right panel shows the 
evolution of each component of energy density for initial conditions <i>o = —7/4 • 5 1//4 , bo = 1/2, 
B = 0. 

V. STABILIZED DILATON 

The vacuum expectation value of the dilaton determines both the gauge and gravita- 
tional coupling constants of the low energy effective theory. Therefore, the dilaton must be 
stabilized at some stage of the evolution for the action to have something to do with the 
reality. In this section, we study the cosmological evolution when the dilaton is stabilized. 
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FIG. 7: Numerical solutions for B ^ and A > 0. When A = 10 -4 , the solutions of four initial 
conditions are given: <j>o = —7/4 • 5 1//4 , bo = 1/2, Bq = for the solid curves, $o = 0, &o = li 
Bq = for the dotted curves, $o = 0, bo = 10, -Bo = for the dashed curves, and <I>o = —^3/2, 
&o = 1, Bq = for the dash-dotted curves. The lower-right panel shows the evolution of each 
component of energy density for initial conditions <I>o = — \/3/2, bo = 1, Bq = 0. 

As for the correct mechanism of dilaton stabilization, the consensus has not been made yet. 
Our goal here is to illustrate an example of the dilaton stabilization and look into the effect 
of it on the dynamics of B field and the cosmological evolution, since the overall features of 
which are insensitive to the detailed mechanism of stabilization. 

Our starting point is the dilaton potential (El). This potential possesses the minimum 
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for A < 0, but the value of the potential at the minimum is negative and need to be set to 
zero by fine-tuning of the constant shift. However, the constant shift of the potential has 
no motivation in the context of string theory. Instead, we introduce a term |m|e -3 * in the 



potential, which can arise from the effect of various form field fluxes in extra dimensions 15]. 
Thus, our potential for the dilaton for B = looks like 

V F ($) = 1 (m|e 3<1> + 2Ae 2 * + m^- 3 *) . (47) 

This potential has a global minimum for any value of A and m F . To have sensible cosmology, 
the potential at the minimum must be zero. For A < 0, this can be done through a fine- 
tuning of the parameters in the potential 

^=l(-i A )"' (48) 

where p? = m F /m B . Then the minimum is located at 

$ F = In (-^ . (49) 

The shape of this fine-tuned potential for A = —0.1 is shown in Figure [HJ Now the equations 
(!35l)-( l39l) are modified accordingly 

13 

$ + (di + a 2 + d 3 )$ = -2p B - p A - -p b - p b + -p F , (50) 
di\ + «i(di + d 2 + d 3 ) = pa + Pb + Pf, (51) 
d 3 + d 3 (di + d 2 + d 3 ) = 2p B + p A + pb + Pf, (52) 

where p F = |/i 2 e _3<1> , while the equation (1541) for b is not changed. 

With the stabilizing potential for the dilaton, the cosmological evolution is completely 
changed. The dilaton and the antisymmetric tensor field rapidly come to the oscillation 



about the potential minimum $ = $f and 6 = 0. Oscillating $ and B fields behave like 
ordinary matter satisfying the equation of state p — The universe becomes isotropic 
and matter dominated. The numerical solutions for the stabilizing dilaton potential (1471) 
with A = —0.1 are plotted in Figure El One can see the oscillation of $ and b, and that 
both di and d 3 approach to 2/3t, indicating the matter domination. Damping of $ and b 
oscillations is due to the expansion of the universe. 
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FIG. 8: (a) The potential V($) for A = —0.1. The minimum of the potential is at $ m = — In 6 
with V($ m ) = 0. (b), (c) and (d) show the numerical solutions of 3>(i), b(t) and dj(t), respectively, 
for A = —0.1 and initial conditions $o = 0, <I?o = and bo = 1. $(t), b(t), and idj(i) approach <E> m , 
0, and 2/3, respectively. 

VI. CONCLUSIONS 

We investigated cosmology of a four-dimensional low energy effective theory arising from 
the NS-NS sector of string theory with a D-brane which contains the dynamical degrees 
of freedom such as the gravity, the dilaton, and the antisymmetric tensor field of second 
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rank, coupling to the gauge field strength living on the brane. The dynamics of the system 
crucially depends on the curvature A and the brane tension tub through which the dilaton 
obtains a potential of the form; Ae 2 * + ^m^e 3 *. Here, the latter becomes the effective mass 
of the antisymmetric tensor field (B- matter). 

In terms of the homogeneous solution in flat spacetime, we first showed how the dilaton 
$(£) and the non- vanishing magnetic component of the tensor field B(t), in one direction 
evolve in time. For positive A, one finds that the dilaton runs away to negative infinity, and 
Bit) reaches a constant value at late time. When A is negative, one can see that the dilaton is 
stabilized at a finite value and then the dilaton (and also B(t)) shows an oscillatory behavior 
around that value, or runs away again depending on the choice of the initial conditions. 
Such a dilaton stabilization, however, produces a negative effective cosmological constant, 
and thus leads to a collapsing universe in string cosmology. 

When B(t) is turned on, the universe undergoes an anisotropic expansion described by 
the Bianchi type-I cosmology. We found the attractor solutions showing the overall features 
of general solutions and confirmed it through numerical analysis. The dilaton <&(£) decreases 
and settles to a logarithmic decrease in time to negative infinity. When the brane tension 
term dominates, the anisotropy is sustained. If there is a positive curvature term, it domi- 
nates finally over the brane tension term as the dilaton rolls down to negative infinity. Then 
the expansion of the universe turns to be isotropic and linear in time. Accordingly, B(t) 
decreases inversely proportional to time. 

There have been various proposals to stabilize the dilaton. In order to study the dynam- 
ics of the B-matter and the stabilized dilaton system, we adopted an option of generating a 
dilaton mass term of the form m 2 r e _3<1> . Then the potential for the dilaton has a global mini- 
mum and the cosmological constant of our Universe can be fine-tuned to a desired value with 
negative A. While the dilaton evolved to a stabilized value, B shows an oscillatory matter- 
like behavior, and the universe expands as in the usual matter-dominated era recovering the 
isotropy. The obtained result is consistent with that of Ref. [3]. 
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